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Derivation of a Fundamental Diagram for Urban Traffic Flow 

Dirk Helbing 

ETH Zurich, UNO Dll, Universitatstr. 41, 8092 Zurich, Switzerland 

Received: date / Revised version: date 

Abstract. Despite the importance of urban traffic flows, there are only a few theoretical approaches to 
determine fundamental relationships between macroscopic traffic variables such as the traffic density, the 
utilization, the average velocity, and the travel time. In the past, empirical measurements have primarily 
been described by fit curves. Here, we derive expected fundamental relationships from a model of traffic 
flows at intersections, which suggest that the recently measured fundamental diagrams for urban flows 
can be systematically understood. In particular, this allows one to derive the average travel time and the 
average vehicle speed as a function of the utilization and/or the average number of delayed vehicles. 



PACS. 89.40.Bb Land transportation 
Kinetic and transport theory of gases 



47.10. ab Conservation laws and constitutive relations - 51.10.+y 



1 Introduction 

The study of urban traffic flows has a long history (see Ref . 
[1] for an overview). For more than a decade now, physi- 
cists have contributed various interesting models, ranging 
from cellular automata [2I3I4J to fluid-dynamic approaches 
|5l6j . Complementary, one should mention, for example, 
Refs. |7l8j as representatives of publications by traffic en- 
gineers, and also continuous microscopic flow models used 
in commercial software tools such as VISSIM. 

One research area in traffic physics is the transition 
from free to congested traffic in urban road networks 
|9ll0j , which started off with the paper by Biham, Middle- 
ton and Levine [11| . Interestingly enough, the spreading 
of congestion seems to share some features with cascading 
failures [T^ . 

In the following, we will focus on the study of funda- 
mental relationships between flow, utilization, and density 
on the one hand and the average velocity or travel time 
on the other hand. Such relationships were already stud- 
ied in the 60ies |i:ill4ll5ll()ll7ll8ll9l2()l21l22l23| . but this 
work primarily took a phenomenological approach. More- 
over, the recorded data in the regime of congested road 
networks scattered enormously. Thus, it was hard to fit 
a curve. Consequently, many different relationships were 
proposed. The probably most wide-spread formula is the 
one published by the Bureau of Public Roads [TB] . Accord- 
ingly, the travel time % of an urban road section i with 
capacity Ci would follow the capacity constraint function 



7^ = 7:" 



(1) 



(see Fig. [5]). Here, Ai is the arrival flow in that road sec- 
tion and T^ the travel time for light traffic conditions. 



while ai w 0.5 and /3i ~ 4 are fit parameters. It is obvious 
that this formula does not diverge when the capacity Ci is 
reached, as one may expect. This is just one of the many 
theoretical inconsistencies of the proposed phenomenolog- 
ical formulas. It is no wonder that the subject of funda- 
mental diagrams for urban traffic has recently been taken 
up again [24,25|26|27|28,29,30,31j . 

The first measurement of an urban fundamental di- 
agram was presented by Godfrey [32]. Ten years later, 
even Robert Herman and nobel prize laureate Ilya Pri- 
gogine addressed traffic flow in cities |33I34| . They tried to 
derive fundamental relationships via a statistical physics 
approach, which however still contained phenomenological 
elements. Their "two-fluid approach" considered moving 
and standing traffic and the mutual intcrdcpendcncies be- 
tween them. 

Recently, the issue of urban gridlock was reanimated 
by Carlos Daganzo [2H]- Together with Geroliminis, he 
presented convincing evidence for a fundamental diagram 
of urban traffic flow [29130] . Before, fundamental diagrams 
were mainly used for the study of freeway systems [55] . 
There, they were very useful to understand capacity ef- 
fects [T] , different traffic states [3T , and the traffic dynam- 
ics, in particularly the propagation of shock fronts [38| . It 
is, therefore, not surprising that people are eager to find 
fundamental relationships for urban traffic as well. 

The most recent progress is an approximation of the 
urban fundamental diagram by Daganzo and Geroliminis, 
which is based on cutting away parts of the flow-density 
plane [31] . It appears possible to construct a relation- 
ship with kinematic wave theory [35] . In contrast to the 
density-based approach by Daganzo et ai, we will pursue 
an alternative, utilization-based approach. This is com- 
mon in queueing theory |39j and transportation planning, 
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where formulas such as the capacity constraint function 
^ are used. 

After discussing elementary relationships for cyclically 
signalized intersections of urban road networks in Sec. [21 
we will start in Sec. [3] with the discussion of undersatu- 
rated traffic conditions and derive a relationship for the 
average travel time as a function of the utilization or the 
number of delayed vehicles. Furthermore, we will deter- 
mine a formula for the average speed. In Sec. IH we will 
extend the analysis to congested road conditions, where 
the intersection capacity is exceeded. Afterwards, in Sec. 
m we will indicate how to deal with oversaturated net- 
works, where the link capacity is insufficient to take up all 
vehicles that would like to enter a road section. Finally, 
Sec.[6]provides a summary and discussion. In particularly, 
we will address issues regarding the transfer of link-based 
fundamental diagrams to urban areas. We also try to con- 
nect the density-based fundamental diagram of Daganzo 
et al. with the utilization-based approach developed here. 



2 Elementary Relationships for Cyclically 
Operated Intersections 

Let us study a single intersection with a periodically oper- 
ated traffic light. We shall have green phases j of duration 
ATj , during which one or several of the traffic streams i 
are served. Pij shall be 1, if traffic stream i is served by 
green phase j, otherwise /3y = 0. The setup time after 
phase j shall require a time period tj . It may be imagined 
to correspond to the time period of the amber light (al- 
though in practice, this has to be corrected for reaction 
times and intersection clearing times) . The sum of setup 
times will be called the "lost service time" 



'los 



3 



while the sum of all green time periods and setup times 
will be called the "cycle time" 



Tcyc = 5] (/IT, +T,)= Tio. + J2 ^T^3 



(3) 



The green times are also sometimes expressed as fractions 
fj > of the cycle time, i.e. 



AT, = f,T, 



eye 



(4) 



with J2i fj < 1- After inserting this into Eq. ^ and re- 
arranging terms, we get 



TcyciifA) 



T, 



los 



1 - E, fj ' 



(5) 



i.e. the cycle time is proportional to the lost service time 

Tlos- 

Assuming average inflows Ai per lane, the number of 
vehicles belonging to traffic stream i, which must be served 
within one cycle time Tcye, amounts to AiTcyc- In order 



to avoid the formation of growing queues (i.e. the onset of 
congestion) , the number of served vehicles per cycle time 
and lane must reach this value. The number of vehicles of 
traffic stream i potentially served during the green phases 
ATj is given by ^jQijPijATj per lane, where Qij de- 
notes the outflow capacity (discharge flow) per lane, when 
stream i is served by phase j. If we want to have a certain 
amount of excess capacity to cope with a variability of the 
inflows, we may demand 



{l + 5^)A,T,, 



j 



^ijPij^-^j / ^ ^ij Pij Jj -^ eye \^) 



with 6i > 0. This linear set of equations may be solved 
for the green time fractions fj. In the following, we will 
assume the case, where not several traffic streams i are 
served in parallel by one and the same green phase, but 
where each phase j serves a single traffic stream i. Then, 



we may assume (3ij = 1 , if j 



and Pi. 



otherwise. 



Furthermore, Qij = Qa = Qi. With this, Eq. ([6]) implies 



where 



A- 
fi{ui, 5i) = (1 + 5i)^ = (1 + 5i)ui , 



A^ 

Q^ 



(7) 



(8) 



is called the utilization of the service or outflow capacity 
Qi (see Fig. [1]). 5^ > is a safety factor to cope with 
variations in the arrival flow (inflow) Ai. According to 
Eq. ([5]) we must have 



0<l-^/, = l-^(l + <5.)^ 



(9) 



Otherwise, we will have growing vehicle queues from one 
signal cycle to the next, corresponding to the congested 
traffic regime. 

Note that, at time t, the number Ni{t) of vehicles per 
lane in the road section reserved for traffic stream i is 
given by the time integral of the arrival flow Ai (t) minus 
the departure flow ji{t)Oi{t): 



N,it)^ dt' A,(i')-7.(i')0.(t') 



(10) 



Here, the starting time io rnust be properly chosen to give 
the correct initial number of vehicles on road section i. 
During the amber and red time periods, we set the per- 
meability 7i(i) = 0, as there is no outflow, while 7i(i) = 1 
during green phases. The departure flow Oi per lane is 
given by the service capacity Qi per lane, as long as there 
is a positive number ANi{t) > of delayed vehicles. If 
Qout represents the characteristic outflow from congested 
traffic per lane into an area of free flow, the overall service 
capacity by all service lanes is given by the minimum of 
the number of lanes U used by vehicle stream i upstream 
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Fig. 1. Schematic illustration of vehicle trajectories for a traffic 
light, which has an amber and red time period of altogether 
(1 — fi)Tcyc and a green time period of fiT^yc- Vehicles move 
forward at the free speed V^P or are stopped in a vehicle queue 
(horizontal lines) , which forms during the amber and red time 
period behind the traffic light (located at the f-axis) . The speed 
of the upstream moving congestion front is given by the arrival 
flow [38] and denoted by Ci < 0. The dissolution speed co < 
of congested traffic is a characteristic constant with jcoj > \ci\ 
[6140) . The average delay time can be determined by averaging 
over the waiting times in the triangular areas. Note that for 
the case of an excess greentime {fi > Ui), vehicles may pass 
the traffic light without any delay. 



the intersection, and the number I^ of lanes downstream 
of it: 



liQi = min{Ii,I-)Qo 



I.e. Qi = min ( 1, y- ) Qout • 



(11) 
When the vehicle queue forming behind a traffic light 
has completely resolved, the outflow from the road section 
used by vehicle stream i drops from Qi to a lower value. 
Then, if the greentime period continues, the outflow Oi{t) 
per lane corresponds to the arrival flow Ai{t — T^^) per lane 
expected at the end of road section i under free flow con- 
ditions [30]. Here, T^ — Li/V^ represents the travel time 
under free flow conditions, which is obtained by division of 
the length Li of the road section reserved for stream i by 
the free speed V.j^ . Considering also the above definition 
of the permeabilities ^i{t) reflecting the time-dependent 
states of the traffic signal, we have 



7.W0.W=7^W 



Qi iiAN,{t)>0, 

Ai{t - 7^") otherwise. 



(12) 



The number ANi [t) of delayed vehicles per lane at time 
t on the road section reserved for vehicle stream i can be 
easily determined as well. In contrast to Eq. ([T0|) . we have 
to subtract the integral of the departure flow from the 
integral of the arrival flow Ai{t — T^) expected at the end 



of road section i under free flow conditions. Altogether, 
the number of delayed vehicles can be calculated as 



AN,{t) = / dt' A,{t' - T^) - j,{t')0,{t') 



(13) 



If the traffic flow is organized as a vehicle platoon and the 
green phase is synchronized with its arrival at the traffic 
light, the number of delayed vehicles is zero. However, if 
the traffic flow Ai is uniform, we find 



AN,{t) =A,-it~to)- / dt' j,{t')0,{t') 



(14) 



to 



Let us assume that to denotes the time when the green 
phase for traffic stream i ended. Then, the next green 
phase for this trafSc stream starts at time ig ^ to + {I — 
fi)Tcyc, as /jTcyc is the green time period and (1 - fi)Tcyc 
amounts to the sum of the amber and red time periods. 
Due to 0^{t) > A,{t) and O^it'^) > A^{t'^), t'^ - to = {1 - 
fi)Tcyc is also the time period after which the maximum 
number AN^'^^ of delayed vehicles is reached. 

In case of a uniform arrival of vehicles at the rate Ai — 
UiQi per lane we have 

ANr^'in^, {/,}) - A,(l - /OTcyc(/0 

= u4.(i-/.)rcyc({/j}). (15) 

Since Qi — Ai is the rate at which this vehicle queue can be 
reduced (considering the further uniform arrival of vehicles 
at rate Ai), it takes a green time period of 



Ti{u^,{Jj}) = 



A,{l-fi)T,y, u,{l-f,) 



Qi Ai 



1 



Jcycd/jl) ■ 



(16) 

until this vehicle queue is again fully resolved, and newly 
arriving vehicles can pass the traffic light without any de- 
lay. 

Finally, let us determine the average delay time of ve- 
hicles. If we have a platoon of vehicles which is served by 
a properly synchronized traffic light, the average delay is 
Tf^ = T™™ = 0. However, if we have a constant arrival 
flow Ai , the average delay 7^**^ of queued vehicles is given 
by the arithmetic mean (7^™^'' + 7^™'")/2 of the maximum 
delay for the first vehicle in the queue behind the traffic 
light, which corresponds to the amber plus red time period 



Tr^'UfA) = i^ - f^)Tcyc{{f,}) , 



(17) 



and the minimum delay 7^™'" = of a vehicle arriving just 
at the time when the queue is fully dissolved. To get the 
average delay, we have to weight this by the percentage 
of delayed vehicles. While the number of vehicles arriv- 
ing during the cycle time Tcyc is AiTcyc, the number of 
undelayed vehicles is given by Ai{ATi —Ti). Considering 
formulas ([4]) and ([T6| . the excess green time is 



^-^ i -^ i — Ji^c^ 



'^i\^ Jij-'-cyc Ji 



1 - Ui 



-T 



(18) 
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Hence, the percentage of delayed vehicles is 



(1 — fi)/{i — Ui) < 1 of vehicles is delayed, together with 
Eqs. ini) and ^ wc find 



(19) 



Altogether, the average delay of all vehicles is expected to 
be 



rriu^AfA) 



(1-/0 rr^'-iUA) (i-/.)'7^cyc({/,}) 



Inserting Eq. ([5]) gives 



1 - u, 



(1-/.)' 



T, 



los 



(1 - ^0 2(1 - E, /,) ' 

while with Eq. ()15|) we obtain 

1 - /, Z\iVf--(u„ {/j}) 



'rriu^,{fA) 



1 - u, 



2uiQ, 



2 

(20) 



(21) 



(22) 



;^ (l-/.)^rcyc({/j-}) 

~ "'^' (1 - «0 2 



(24) 



Therefore, the average delay time is proportional to the 
maximum queue length AN^^^, but the prefactor de- 
pends on /i = (1 + Si)ui (or the safety factor 5i, re- 
spectively). In case of no excess green time {6i — and 
fi = Ui), we just have 



This remarkably simple relationship is known in queuing 
theory as Little's Law pS], which holds for time-averaged 
variables even in the case of non-uniform arrivals, if the 
system behaves stable (i.e. the queue length is not sys- 
tematically growing or shrinking) . It also allows one to es- 
tablish a direct relationship between the "average vehicle 
density" 

AN"" 

pT - -r- (25) 

on the road section used by stream i and the utilization 
Ui = Ai/Qi, namely 



w»Q,^" 



(26) 



rr[{u,}) = {i-u.) 



Teyc({%}) _ ANr-{{u,}) 



2A, 



2uiQi 



(23) 



(see Fig. O. 




Note that the average density pf has the same dependen- 
cies on other variables as Z\ A^f ^ or Tf-^ , and an additional 
dependency on Li, i.e. the more natural quantity to use is 
the queue length Z\Af^. 



2.1 Efficiency of Traffic Operation 

In reality, the average delay time will depend on the time- 
dependence of the inflow Ai (t) , and on how well the traffic 
light is coordinated with the arrival of vehicle platoons. In 
particular, this implies a dependence on the signal offsets. 
In the best case, the average delay is zero, but in the worst 
case, it may also be larger than 



^-ut .. , _ (1 - -»,)rio, 

2 ^"'^^"^^^"2(1-E,".) 



(27) 



see Eq. (|20|) with fi = m. We may, therefore, introduce 
an efficiency coefficient e^ by the definition 



Tr{{u,},ei) = {l-e,) 



(1 - Ui)T\os 



Fig. 2. Schematic illustration of vehicle trajectories and signal 
operation in the case fi — Ui, where there are no excess green 
times so that the traffic light is turned red as soon as the vehicle 
queue has been fully resolved. 



Similarly to the average travel time, we may deter- 
mine the average queue length. As the average num- 
ber of delayed vehicles is [ANf^'^^ -f 0)/2 and a fraction 



2(1 -E,".) 

or, considering Eq. (j24p . equivalently by 

(1 - Ui)T\o. 



ANr{{uj},ei) = {l-ei)u,Q, 



2(1 -E,%) 



(28) 



(29) 



For ei — 1, the traffic light is perfectly synchronized with 
platoons in traffic stream i, i.e. vehicles are served with- 
out any delay, while for e^ = 0, the delay corresponds to 
uniform arrivals of vehicles, when no excess green time is 
given. If the traffic light is not well synchronized with the 
arrival of vehicle platoons, we may even have e^ < 0. It 
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also makes sense to define an efficiencies not only for the 
traffic phases, but also for the operation of a traffic light 
(i.e. the full cycle). This can be done by averaging over 
all efficiencies e^ (and potentially weighting them by the 
number UiQiTcyc of vehicles arriving in one cycle. There- 
fore, it makes sense to define the intersection efficiency 
as ^ 



Y^i^iQi 



(30) 



Note that, particularly in cases of pulsed rather than uni- 
form arrivals of vehicles, the efficiencies ti depend on the 
cycle time Tcyc and, therefore, also on the utilizations Ui. 
Increasing the efficiency e^ for one traffic stream i will of- 
ten (but not generally) reduce the efficiency Cj of another 
traffic stream j, which poses a great challenge to traffic 
optimization. 

The exact value of the efficiency e; depends on many 
details such as the time-dependence of the arrival flow 
Ai{t) and its average value Ai, the length Li of the road 
section, and the signal control scheme (fixed cycle time 
or not, adaptive green phases or not, signal offsets, etc.). 
These data and the exact signal settings are often not fully 
available and, therefore, it is reasonable to consider e^ as 
fit parameters rather than deriving complicated formulas 
for them. Nevertheless, we will demonstrate the general 
dependence on the utilization Ui in the following. 

For this, we will study the case of excess green times 
{6i > 0), which are usually chosen to cope with the 
stochasticity of vehicle arrivals, i.e. the fact that the num- 
ber of vehicles arriving during one cycle time is usually 
fluctuating. The choice Si > 0, i.e. fi > Ui, also implies 



AT ^ ^ 



T, 



eye 



ATi 
T, 



''-Q^■ 



(31) 



eye 



This reflects the well-known observation that the average 
departure flow Ai = UiQi usually does not reach the value 
given by the green time fraction fi times the saturation 
flow Qi. 

The efficiency ti related with a value 5i > Q may be 
derived from Eqs. ([21]) and (HH). We obtain 



(32) 



1 - e. = 



{i-u,Y (i-E,/.) 



where fi{ui,Si) = {1 + Si)ui according to Eq. ([7|). The 
efficiency e; is usually smaller than in the case, where the 
traffic light is turned red as soon as a vehicle queue has 
been dissolved (for exceptions see Ref. 41 ). Then, Ci < 
for 5i > 0. Formula (|5^ also allows one to treat the 
case where the green time fractions fi and the cycle time 
Tcyc are not adapted to the respective traffic situation, 
but where a fixed cycle time Tcyc = Tc^yc and fixed green 
time fractions /? are implemented. This corresponds to 
constant green times 



/fTe"ye({/n) 



/°rios 



(33) 



In case of uniform vehicle arrivals, we just have to insert 
the corresponding value fi — fi into Eq. ([311) to obtain 
ei. In the case of non- uniform arrivals, e^ can be under- 
stood as fit parameter of our model, which allows us to 
adjust our formulas to empirical data and to quantify the 
efficiency of traffic light operation. In this way, we can 
also absorb effects of stochastic vehicle arrivals into the 
efficiency coefficients ei, which simplifies our treatment a 
lot. 



3 Fundamental Relationships for 
Undersaturated Traffic 

The travel time is generally given by the sum of the free 
travel time 7^" = Li/V^ and the average delay time 7^**^, 
where Li denotes the length of the road section used by 
vehicle stream i and V^ the free speed (or speed limit). 
With Eq. 1241), we get 



(34) 
Inserting Eq. ([29l) . we can express the travel time solely 
in terms of the utilization Ui, and we have 

^(w.^j-^ + d-o J,^:^^;-, ^ (35) 

The formula ([35l) constitutes a fundamental relationship 
between the average travel time 7^*'^ on the capacity uti- 
lization Ui under the assumptions made (mainly cyclical 
operation with certain efficiencies ei). Of course, one still 
needs to specify the factor (1 — e^). In case of constant 
arrival rates Ai, this factor is given by Eq. ([32| . which 
finally results in 

^(«.,ftt)^^ + ^5jr^^ (36) 

After insertion of Eq. (O , we get 

T-.r XSXW ^^ I [1- (l+(5,)ui]^Tios 

^({%}, m - ^ + (i-.02[i-E,(i+^.K] • 

(37) 
Sometimes, it is desireable to express the fundamental 
relationships in terms of the density rather than the util- 
ity. Inserting Eq. ^ into Tt''{ui,ei) = %{ui, ei)-L,/V^, 
and this into Eq. ((251), we obtain the equation 

p, (K},e„T.) = -^(l-.,)2(Y^^-^, (38) 

which can be numerically inverted to give the utilization 
Ui as a function of the scaled densities pYLj/{l — ej). The 
calculations are simpler in case of a fixed cycle time T?,^ 



and an uniform arrival of vehicles. By inserting Eq. ([2(1 
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into dMl) [and with /, = /f, Teyc = T^^yc see Eq. §^], we 
get 



= UiQ 
which finally yields 



— w^^) — ' ^ ^ 



This can be inserted into Eq. (|26p to give 

Pi ^i 



(40) 



rr{pTUAm) = 



u,(pri.,{/,°})Q 



^^ + (1-/°)' ^^ ' ■ (41) 



<9 



3.1 Transition to Congested Traffic 

The utilizations Ui increase proportionally to the arrival 
flows Ai^ i.e. they go up during the rush hour. Eventually, 



E/j = E(i + ^^K^i 



which means that the intersection capacity is reached. 
Sooner or later, there will be no excess capacities any- 
more, which implies (5j — > and ji — > Ui. In this case, 
we do not have any finite time periods ATi — Ti, during 
which there are no delayed vehicles and where the depar- 
ture fiow Oi{t) agrees with the arrival flow Ai. Therefore, 
Oi{t) = qi{t)Qi according to Eq. (fT2|) . and certain rela- 
tionships simplify. For example, the utilization is given 
as integral of the permeability over one cycle time T^yc, 
divided by the cycle time itself: 



t.o+T„ 



1 



Tr, 



dt'j^it') 



(43) 



Moreover, in the case of constant in- and outflows, i.e. lin- 
ear increase and decrease of the queue length, the average 
number ANf^ of delayed vehicles is just given by half of 
the maximum number of delayed vehiclesu 



ANf 



ANf 



(44) 



^ Note that the formulas derived in this paper are exact un- 
der the assumption of continuous flows. The fact that vehicle 
flows consist of discrete vehicles implies deviations from our 
formulas of upto 1 vehicle, which slightly affects the average 
travel times as well. In Fig. [2j for example, the number of 
vehicles in the queue is iV™'"' - 1 rather than iV™''''. There- 
fore, the related maximum delay time is (1 — Ui)Tcyc{N'^'^^ — 
1)/N^'", which reduces the average delay time T^" by (1 - 



As a consequence, we have 






(45) 



3 "■!) 



The average speed V^-'^ of traffic stream i is often deter- 
mined by dividing the length Li of the road section re- 
served for it by the average travel time % = Tl' + T^^ , 
which gives 



vr{{u,}) - 



1 , {l-Ui)T,, 



%{{u,}) IVP 2L,(1-E 



j^]) 



(46) 

and can be generalized with Eq. (pS)) to cases with an 
efficiencies ti ^ Q: 

(47) 
Note, however, that the above formulas for the average 
speed are implicitly based on a harmonic rather than an 
arithmetic average. When correcting for this, Eq. (|1^ . for 
example, becomes 



(42) Vr{{u,}) - 



J-^i'^i'^i 



^^max(|y^.|) 



In 



zAiV-- ({«,}) 



i^i -^i 



As is shown in Appendix |^ this has a similar Taylor ap- 
proximation as the harnronic average (|47p . The latter is 
therefore reasonable to use, and it is simpler to calculate. 
As expected from queuing theory, the average travel 
time (|45p diverges, when the sum of utilizations reaches 
the intersection capacity, i.e. ^ 



J "J 



1. In this practi- 
cally relevant case, the traffic light would not switch any- 
more, which would frustrate drivers. For this reason, the 
cycle time is limited to a finite value 



Tio 



ymax/ r Ul.\ ^ ^^os 



(49) 



3 J 



where typically m^ < Uj. This implies that the sum of 



utilizations must fulfill 



E-.^E"" = i 



T, 



los 



'T'max. 
eye 



(50) 



As soon as this condition is violated, we will have an in- 
crease of the average number of delayed vehicles in tinre, 
which characterizes the congested regime discussed in the 
next section. 



4 Fundamental Relationships for Congested 
Traffic Conditions 

In the congested regime, the number of delayed vehi- 
cles does not reach zero anymore, and platoons cannot 
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be served without delay. Vehicles will usually have to 
wait several cycle times until they can finally pass the 
traffic light. This increases the average delay time enor- 
mously. It also implies that there are no excess green 
times, which means 6i = 0. Consequently, we can also as- 
sume Oi{t) = Qi, as long as the outflow from road sections 
during the green phase is not (yet) obstructed (otherwise 
see Sec. [5]). Formula (IT51) applies again and implies for 
time-independent arrival flows Ai 



tio+kT^yT 



AN,{t,o + kT^:n ^ ^N^iUo 



dt' 



A^~j^{t')Q, 



AN,(Uo) 



tiO 



where 



However, this turns out to be not correct, as some of this 
time is actually not lost, but used to move forward (see 
Fig. [3]). After this flrst simple analysis, we need to look at 
the problem more carefully, and perform the calculation 
without the use of Eq. ([M]) . as Little's Law is not applica- 
ble in case of non-stable, systematically growing queues. 

Again, we apply the fact that the number of vehicles 
arriving in one cycle is AiT^^^, while the number of ve- 
hicles served during one green time period of duration 
u,°T™j^'' is given by QjM^T™^^''. The difference of both num- 
bers is added to the growing vehicle queue. The quotient 
of the number kAiT"^^^ of vehicles arriving in k cycles 

and the number Qiu'^T^^^ of vehicles served during one 
green time period, when rounded down, corresponds to the 
number rig of additional stops needed by newly arriving 
vehicles. Therefore, with Ai — uiQi we have 



Uo+kT" 






dt' 7,(t') < u. 



A. 



(52) «s(ui,({u°}),/e) = 



tiO 



is the green time fraction of the cycle time T^y^^ re- 
served for vehicle stream i. Therefore, the number of de- 
layed vehicles grows by an amount {Ai — u^Qi)T^l^ ~ 

{ui — u^)QikT^y^^ in each cycle time T^y^^- The minimum 
number during one cycle is 



Z\iVr"(M„ k) = AN,{t,o) + [u^ ~ u1)QikT, 



max 
eye : 



(53) 



and considering Eq. (jlSp , the maximum number of delayed 
vehicles is 



AN^^iu,, k) = Z\iVr"(w„ k)^u,{l-u1)Q,T, 



max 

eye ■ 



(54) 



Because of ANf^ = (ziA^f"" + AN^'^'')/2 and A, ^ u^C 
the average number of delayed vehicles is 



Z\iVr(u„{w°},fc)= [{u.~u°)k 



+ ^^'K^]QrT^--{{u]}). (55) 



Obviously, the average density 

7..(i-wO)\g.r-f(^o}) 



u 



(56) 






(HI) 



Q,,07^max 
I "i -^ eye 



({«?}) 



Uik 



(58) 



Here, \_x\ represents the floor function rounding to lower 
integers (i.e. cutting the digits after the period, if a; > 0). 
Assuming that i^o is the time at which congestion sets in, 
Eq. (|58p can be generalized to continuous time t, allowing 
us to estimate the number of additional stops of a vehicle 
arriving at time t: 



ns{u^,{u°A,t) 



A,it~tio) 



■ ii^T'niax 



(Kl) 



({^°}) 



.OT^niax 
■''i eye 



(59) 

As Fig. [3] shows, the delay time by each of the Us ad- 
ditional stops is (1 — u'^)T^y^'^. Moreover, we can see that 
the triangular part of the vehicle queue in the space-time 
plot gives a further contribution to the delay time of ve- 
hicles. Applying Eq. pO|) with fi — m = u'-, the average 
time delay in this triangular part is (1 — u'^)T^^^/2, i.e. 
the arithmetic average between zero and the sum of the 
red and amber time period (amounting to (1 — u'^)T^y^^). 
In summary, the delay time of a newly arriving vehicle at 
time t (when averaging over the triangular part for the 
sake of simplicity), is 



Triu,,t) = {i-u'^)^^ + 



T^inax 
eye 



Ui{t ~ tio) 



^.OT^max 

"i -^ eye 



max 

eye 



Uiit-tio) 



■^i eye 



{i-u'1)t^;T- (60) 



is obtained by dividing the previous formula by Li. 

It seems logical to use Eq. ([55)1 to determine the aver- 
age delay time as 



7:-(^„{<},fc) = 



Z\7Vf"(u„{wO},fc) 

A. 



1 - - I kT^yTiiu'j}) 
(1 - 7.0)T-r({7.0}) 



(57) 



Accordingly, the average travel time does not only grow 
with time t (or the number k of cycles passed), it also 
grows stepwise due to the fioor function \x\ . 

When we also average formula (|60p over its steps, we 
obtain the approximate relationship 



Triu..,t) 



Uiit - Uo) 



^.OT^max 
"'i eye 



(i-"?)T, 



max 

eye 



Uj(i- tjo) 



(1 



(61) 
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Fig. 3. Schematic illustration of growing vehicle queues in the 
congested regime. Since not all vehicles can be served within 
one green time u^T^f^, vehicles are forced to make additional 
stops, which correspond to the rhomboidal parts of the queue in 
the above time-space diagram. The triangular part, in contrast, 
represents newly arriving vehicles. 



where it is important to consider that the floor function 
[a;J is shifted by 0.5 with respect to the function round(a;), 
which rounds to the closest integer: round(a;) — [a; + 0.5j. 
Moreover, when averaging over the last (fc + 1st) cycle we 
get 

/ 1 \ 1 - u° 

rr{u., k) ^u,[k + - j ^^T-r . (62) 

It is also possible to replace the dependence on the num- 
ber k of cycles by a dependence on the average density 
of delayed vehicles by applying Eq. ([SS)) . In this way, we 
obtain 



^.\- 



ANr{u,,{u^},k) u^,{l-Ui) 



{m - ^i°)Q.r--({«0}) 2iu, - u^) 



(63) 



Therefore, while in Sec. [3] we could express the average 
travel time and the average velocity either in dependence 
of the average density pf or the utilization Ui alone, we 
now have a dependence on both quantities. 

Finally note that Eqs. ([SO)) to (|62p may be generalized 
to the case where the arrival rate of vehicles is not time- 
independent. This changes the triangular part of Fig.[3l In 
order to reflect this, the corresponding contribution (1 — 
ii°)T™^'^/2 may, again, be multiplied with a prefactor (1 — 
Ci), which defines an efficiency e^. While ei — corresponds 
to the previously discussed case of an uniform arrival of 
vehicles, e^ = 1 reflects the case where a densely packed 



platoon of vehicles arrives at the moment when the last 
vehicle in the queue has started to move forward. 



5 Fundamental Relationships for 
Oversaturated Traffic Conditions 

We have seen that, under congested conditions, the num- 
ber of delayed vehicles is growing on average. Hence, the 
vehicle queue will eventually fill the road section reserved 
for vehicle stream i completely. Its maximum storage ca- 
pacity per lane for delayed vehicles is 



ANriu) = i,pi 



jam 



(64) 



where p'^^"^ denotes the maximum density of vehicles per 
lane. The road section becomes completely congested at 
the time 

t == t,o + fcr"r + ^^ (65) 



when AN^'''{ui,k) + AiAt according to Eq. ([53]) reaches 
the value AN^'^"^, which implies 



At 



ANl' 



ANf^'^u^.k) 



where 



fc = 



A, 



ANf 



{Ui 



,0\T^max 
z / eye 



(66) 



(67) 



see Eq. ([55)) . The number n^ + 1 of stops is given by Eq. 
([59)) . See Fig. d for an illustration. 

Complete congestion causes spillover effects and ob- 
structs the arrival of upstream vehicles, even though the 
green phase would, in principle, allow them to depart. 
When these spill-over effects set in, we are in the over- 
saturated regime, and only a certain fraction Ui of the 
green phase of duration u^T^f^ can be used, where < 
<Ji < 1. Note that Ui decreases in time as the number of 
blocked subsequent road sections grows. This may even- 
tually lead to gridlock in a large area of the urban traffic 
system. 

In the case a^ < 1, fi — u^ must essentially be re- 
placed by fi = aiu^ in the fundamental relationships for 
congested traffic, and the lost service time becomes 



Tios = E h + (1 - ^■')"^^< 



OT^inax 
eye 



(68) 



That is, the reduced service time may be imagined like an 
extension of the amber time periods. Therefore, it would 
make sense to reduce the cycle time to a value Tcyc < T^yf^ 
in the oversaturated regime. 

Note, however, that the travel times on the road sec- 
tion reserved for stream i are not growing anymore, be- 
cause the road section is limited to a length Li. This allows 
us to determine the corresponding travel time on the link 
as follows: The number of vehicles served per cycle time is 
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Fig. 4. Schematic illustration of the service of vehicle queues, 
when the road section is fully congested. The lower horizon- 
tal line indicates the location of the upstream end of the road 
section. It can be seen that vehicles are stopped several times, 
and that new vehicles can only enter when some vehicles have 
been served by the traffic light located at the i-axis, and the 
space freed up by this has reached the end of the vehicle queue. 
For similar considerations see Ref. [35) . Note that the charac- 
teristic speed Co of the shock fronts (which corresponds to the 
slope of the congested flow-density relationship for a road sec- 
tion without the consideration of traffic lights) is different from 
the slope of the urban fundamental diagram, because of the ef- 
fect of signal offsets and delays ^3JL. Therefore, co should be 
understood as fit parameter, here. 



aiU^T^'^^Qi. For this reason, the average travel time can 



be estimated as 



li{Ui,ai,Li) = TTTTr——!^, 



Uf^t 



^ ,,,07imax eye q^ i 



and the average delay time TJ^'^ = 7^ — T^'' as 



(69) 



Tr{u^,a,,U) 



T n*' 



L, 



a.u'iQ, V^ 



(70) 



Note that these values are now independent of both, the 
utilization and the average density, as soon as the latter 
assumes the value pf — p^/™, corresponding to a fully 
congested road section. 



5.1 Transition from Oversatu rated to Undersaturated 
Traffic Conditions 

If the arrival flow Ai after the rush hour drops below 
the value of am^T^f^, the vehicle queue will eventually 
shrink, and the road section used by vehicle stream i en- 
ters from the oversaturated into the congested regime. The 
formulas for the evolution of the number of delayed vehi- 
cles are analogous to Eqs. ([55]) and ([5^ . The queue length 
starts with ANl'^"^ and is reduced by {ui—aiU^)QikTcyc < 
in each cycle of length Tcyc. Counting the number of cy- 
cles since the re-entering into the congested regime by k' , 
we have 

AN^'{u^, k') = ANj^"' + (u, - a,u°)Q,k'T,y, , (71) 



and considering Eq. ()15p , the maximum number of delayed 
vehicles is 

Z\7Vf-'^(u„ A:') = ANf^iu,, k') -f u,{l - <J,u°)Q,T,y, . 

(72) 
As soon as AN™^'^{ui,k') reaches zero, the road section 
used by vehicle stream i enters the undersaturated regime. 
Before, the number of stops of vehicles joining the end of 
the vehicle queue are expected to experience an number 
Us + 1 oi stops with 



ris(fc') = 



(73) 



compare Eq. 



6 Summary and Outlook 

Based on a few elementary assumptions, we were able to 
derive fundamental relationships for the average travel 
time T^^ and average velocity V]^^ . These relationships 
are functions of the utilization Ui of the service capacity 
of a cyclically signalized intersection and/or the average 
number ANf''^ of delayed vehicles (or the average density 
p^^ of vehicles in the road section of length L, reserved 
for traffic stream i). We found different formulas, (1) for 
the undersaturated regime, (2) for the congested regime, 
and (3) for the oversaturated regime. While we also dis- 
cussed situations, where fixed cycle times T°y^ are applied, 
we primarily focussed on situations, where the cycle time 
is adjusted to the utilization Ui (in the undersaturated 
regime) and to the effectively usable green time fraction 
ffi (in the oversaturated regime). Our results are summa- 
rized in Fig. [5]), where also a comparison with the capacity 
restraint function ^ is made. 

The formulas for the non-congested regime can be ei- 
ther expressed as non-trivial functions of the utilization 
Ui or the average queue length ANf^ (or the average den- 
sity p^^^)- They contain a fit parameter e^, which reflects 
effects of variations in the arrival ffow and relates to the 
efficiency of traffic signal operation in terms of synchro- 
nizing with vehicle platoons. In the best case, delay times 
are zero, which shows the great optimization potential for 
traffic control in this regime. 
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0.4 0.6 

Utilization Uj 

Fig. 5. Schematic illustration of the capacity restraint func- 
tion ([T]) that the Bureau of Public Roads recommends to use 
|16) (solid line), together with analytical results of this paper 
(dashed lines). Travel time % is measured in units of T^ . The 
green long-dashed line shows that the travel time diverges at 
Ui ~ 0.45, if two green phases during one cycle and identical 
utilizations ui = U2 of both associated road sections are as- 
sumed, furthermore, if the parameters are set to Ties = 0.17^" 
and S = 0.1. (A generalization to traffic operation with more 
than two green phases is easily possible.) When the cycle time 
is limited to a finite value 7^™^" to avoid infinite delay times in 
one of the vehicle queues, one will have growing vehicle queues 
and increasing travel times in both road sections. Therefore, 
the travel time can assume any value above the lower dashed 
horizontal line. The link travel time is only limited by the 
above dashed horizontal line, which corresponds to the situ- 
ation where the vehicle queue fills the road section completely. 
Note that the capacity restraint function (solid line) averages 
over all travel time measurements for a given utilization Ui. In 
the right part of the diagram this concerns measurements that 
depend on the duration of congestion and scatter between both 
horizontal dashed lines. Here, the curve is shown for ai — 0.5, 
f3i = 4, and A,/a = u,/0.45. 



In the congested regime, the number of delayed vehicles 
grows in time, and the majority of vehicles is stopped sev- 
eral times by the same traffic light. Therefore, the average 
travel time does not only depend on the utilization Ui, but 
also on the average vehicle queue ANf^ (or the average 
density pf). Although the traffic light control can still im- 
prove the average travel times by synchronizing with the 
arrival of vehicles, the related efficiency effect is rather 
limited. 

In the oversaturated regime, the storage capacity of the 
road section is fully occupied by delayed vehicles, which 
obstructs the arrival flow. Therefore, the average travel 
time of a road section reaches a constant maximum value. 
Nevertheless, the travel time of vehicles increases further 
in time due to spillover effects, which trigger the spread- 
ing of congestion to upstream road sections. The actually 
usable fraction of the green time period is described by 
a parameter Ui. Synchronization can still reach some im- 



provements. The most favorable control is oriented at a 
fluent upstream propagation of the little remaining free 
space (the difference between ANf'''^ and Z\iV™'"). That 
is, rather than minimizing the delay of downstream mov- 
ing vehicle platoons, one should now minimize the delay 
in filling upstream moving gaps [42j . Furthermore, note 
that the free travel time T^ = L-JV^ and the delay time 
in the oversaturated regime are proportional to the length 
Li of the road section, while the delay times in the under- 
saturated and congested regimes are independent of Li . 

Based on the above results, it is obvious that it cannot 
be very successful to describe the travel time of a link by 
a capacity restraint function which depends on Ai/Ci = 
Ui/u^ only. This basically means an averaging over data 
that, in principle, are also dependent on the average queue 
length AN'f (or on the time passed since the onset of 
congestion). It is, therefore, no wonder that empirical data 
of travel times as a function of the utilization scatter so 
enormously in the congested and oversaturated regimes 
(see e.g. Ref. [27]), that a fitting of the data to functional 
dependencies of any kind does not make much sense. 

This has serious implications for transport modeling, 
as capacity restraint functions such as formula ([1]) are used 
for modeling route choice and, hence, for traffic assign- 
ment. Based on the necessary revision of this formula and 
comparable ones, all traffic scenarios based on these for- 
mulas should be critically questioned. Given the computer 
power of today, it would not constitute a problem to per- 
form a dynamic traffic assignment and routing based on 
the more differentiated formulas presented in this paper. 



6.1 Transferring the Link-Based Urban Fundamental 
Diagrams to an Area-Based One 

We may finally ask ourselves, whether the above formu- 
las would also allow one to make predictions about the 
average travel times and speeds for a whole area of an ur- 
ban traffic network, rather than for single road sections 
("links") only. This would correspond to averaging over 
the link-based fundamental diagrams of that area. For 
the sake of simplicity, let us assume for a moment that 
the parameters of all links would be the same, and de- 
rive a velocity-density diagram from the relationship ()85|) 
between average vehicle speed V^"^ and the capacity uti- 
lization Ui. Taking into account V,^ = Li/T^ and dropping 
the index i, we can write 

^.av(^)^ ^°l^{l+[l-^/(^.)]^eyc/r"} ^ ^^,f{u)-u 



(1 - u)T,y,/TO 
with f{u) = (1 + d)u. The cycle time 



1-u 



T, 



Tlos 



eye 



1 - sfiu) 



(74) 



(75) 



follows from Eq. ([5]) , assuming s signal phases with fi — f 
for simplicity. Note that the previous formulas for pf^ de- 
note the average density of delayed vehicles, while the av- 
erage density of all vehicles (i.e. delayed and freely moving 
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ones) is given by 



p{u) = 



N{u) _ N{u)/T(u) _ A{u) 



L 



L/T{u) ^"(u) V^{u) 



(76) 



where N = AT = uQT denotes the average number of 
vehicles on a road section of length L. Plotting V^^{u) 
over p{u) finally gives a speed-density relationship V'^^{p). 

We have now to address the question of what happens, 
if we average over the different road sections of an urban 
area. Considering the heterogeneity of the link lengths i^, 
efficiencies e^, utihzations u^, and densities pf^ , one could 
think that the spread in the data would be enormous. 
However, a considerable amount of smoothing results from 
the fact that in- and outflows of links within the studied 
urban area cancel out each other, and it does not mat- 
ter whether a vehicle is delayed in a particular link, or 
in the previous or subsequent one. Therefore, the result- 
ing fundamental diagrams for urban areas are surprisingly 
smooth 29t30j. The details of the curves, however, are ex- 
pected to depend not only on the average density, but also 
on the density distribution, the signal operation schemes, 
and potentially other factors as well. 

When averaging over different links, we have to study 
the effects of the averaging procedure on the density and 
the speed. The density just averages linearly: Thanks to 
Little's law [35], the formula ([7S|) can also be applied to 
an urban area, as long as the average number of vehicles 
in it is stationary. However, as the link-based fundamental 
diagram between the flow Q{p) — pV^^{p) and the density 
p is convex, evaluating the flow at some average density 
overestimates the average flowo This also implies that, 
when the relationship V^'^{p) is transferred from single 
links to urban areas, the average speed is overestimated 
for a given density, as Ref. [3T] has shown. 

Despite this expected deviation in heterogeneous and 
inhomogeneously used road networks, and despite the 
many other simplifications, the curve V'^'^{p) fits empir- 
ical data of the speed-density relation in an urban area 
quite well. Figure [6] displays empirical data obtained for 
the center of Yokohama [521 together with a fit of the the- 
oretical speed-density relationship y^^(p), where only the 



^ When averaging over speed values, they have to be 
weighted by the number of vehicles concerned, i.e. by the den- 
sity. This comes down to determining an arithmetic average 
of the flow values and dividing the result by the arithmetic 
average of the related densities. Hence, an overestimation of 
the average flow also implies an overestimation of the aver- 
age velocity. However, considering the curvature of Q{p) and 
knowing the variability of the density p allows one to estimate 
correction terms. 

^ This corresponds to 10% extra green time, an average dis- 
tance between successive traffic lights of roughly 100 meters 
(depending on Tios), and an average of 3 traflic phases (which 
appears plausible, considering that there are many one-way 
roads, which need less than 4 phases in one cycle time). Con- 
sequently, all parameters are quite reasonable (see also Ref. 
[31)). Note that efi'ects of oversaturation did not have to be 
considered in Fig. (6] This is, in fact, consistent with pictures 
from Google Earth. 




0.04 0.06 
Density (1/m) 



Fig. 6. Fundamental velocity-density relationship for a central 
area of Yokohama. Small circles correspond to empirical data 
by Kuwahara as evaluated by Daganzo and Geroliminis [31) . 
The fit curve corresponds to the theoretically derived equa- 
tions (|74p to (|76p . where the outflow (discharge flow) Q = 
1800 veh. /hour/lane and the free speed 1/° — 50 km/h have 
been fixed. The only fit parameters were 5 — 0.1, Tios/To = 1-4, 
and s = 3.^ 



three parameters Vq, (5, and V* = L/T\os were adjusted. 
Surprisingly, the effects of network interactions could be 
sufficiently well represented by a single parameter S, which 
relates to the efficiency e of road sections according to Eq. 
P2p . This approximation seems to work in situations close 
enough to a statistical equilibrium (when the number of 
vehicles in the urban area does not change too quickly) . 

In contrast, for an understanding of the spreading dy- 
namics of congestion patterns, we expect that one must 
study the interaction between the flow dynamics and the 
network structure. This difficult subject goes beyond the 
scope of this paper and beyond what is doable at the mo- 
ment, but it will be interesting to address it in the future. 
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A Determination of Average Travel Times 
and Velocities 

Let /(x) be a function and w{x) a weight function. Then, 
the average of the function between x ~ xq and a; = a:i is 
defined as 



over wc/a from xi = AN^"""" to 2:2 = (1 - Ui)A.,T'^y'= = 
{l — Ui)uiQiT'^y'^, where the specifications of a, 6, and c are 
unchanged. Taking into account V^^ = Li/T^ , this gives 



/ dx' w{x')}{x') 

xo 

xi 

J dx' w{x') 



(77) 



where 



'jL,u,Q, In |1 + AN^^^/{u,Q,T^)\ + Z 



(83) 



In case of uniform arrivals of vehicles, we have a functional 
relationship of the form f{x) = a + bx for the travel time, 
and the weigth function is constant, i.e. w{x) = w. Here, 
a = 7;", X = AN,, and b = 1/A, = 1/(m,Q,). With xo = 
and xi = AN™'^^, the formula for the average travel time 
becomes 



Z = ti;y,0[(l - u,)u,Q,T''y'' - ANr""] ■ (84) 

Considering Eq. (fTS)) . we get 



V 

' 



u 



(1 - u.;)T, 



eye 



In 1 + (1 - /,) 



Tr, 



rp 



V, 



h ~ Ui 



1 — Uj 



In second-order Taylor approximation, this results in 



(85) 



[{axi + bxi^/2) - {axp + bxp^ /2)] 
w- {xi- Xo) 



= a + b 



Xi + Xq 






' 1 



1 - U,j 



(1 ~ /^)Te. 

2T° 



fi -Mj 



1 - Ui 



(86) 



(78) 

where we have used (xi^ — xq^) — (2^1 — a;o)(2;i + xq). 
Inserting the above parameters, we obtain the previously 
derived result 

7:^7-" + ^^!^. (79) 

2uiQi 

When determining the average velocity Vf''^ , the function 
to average over is of the form f{x) = c/(a + bx), where 
c = Li and the other parameters are as defined before. We 
use the relationship 



which can also be derived from Eq. ([5^ . considering Eq. 
(fT5| and the percentage of delayed vehicles, which is given 
by Eq. P^ . The same result follows from 



vr ^ 



together with Eq. ([2T 



T ■ / 7"^v 



(87) 



dx' 



a + bx' 



— { In |a + bxi I — In |a + bx^ 



wc , 



a + bxi 



a + bxo 



(80) 



Dividing this again by the normalization factor w-{xi —xq) 
and inserting the above parameters finally gives 



' 1 






In 



ANY 



u,Q2rp 



L, / ZiiV"'^'^ 



(81) 

where we have used ln(l + x) < x — x^ 12. This formula 
corrects the naive formula 




]/av ~ _i = 



luiQi 



^^max 



(82) 

where we have used 1/(1 + x) « 1 — x. Therefore, the 
above Taylor approximations of both formulas agree, but 
higher-order approximations would differ. The formulas 
in the main part of the paper result for Nf^ = Nf^"^^ /2, 
which corresponds to the case i5i = (i.e. fi — Ui). 

Generalizing the above approach to the case Si > 0, 
we must split up the integrals into one over wc/{a + bx') 
extending from xq = to xi = AN^^'^^ and another one 



